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The first order electroweak phase transition in the Standard Model turns into a regular cross-over 
at a critical Higgs mass mu.c ~ 75 GeV. At the critical point the transition is of second order. 
We make a detailed investigation of the critical properties of the electroweak theory at the critical 
point, and we find that the transition falls into the 3d Ising universality class. The continuum 
limit extrapolation of the critical Higgs mass is iriH,c = 72(2) GeV, which implies that there is no 
electroweak phase transition in the Standard Model. 
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I. INTRODUCTION 

During the last several years, an impressive amount of 
quantitative knowledge about the Standard Model (SM) 
finite temperature electroweak phase transition has been 
obtained through lattice Monte Carlo simulations. At 
small Higgs masses {run) the transition is of the first 
order. The transition becomes weaker (decreasing latent 
heat, surface tension) when the Higgs mass increases, and 
it has been found to turn into a regular cross-over when 
niH > 75 GeV • At the endpoint of the first order 
transition line a second order phase transition appears. 
At this point the macroscopic behaviour of the system is 
determined by the universal properties of the endpoint. 
While the location of the endpoint and the mass spec- 
trum near it have been studied before, the critical prop- 
erties of the endpoint itself have not been resolved so far. 
We shall show that the universality class of the SM end- 
point is of the 3d Ising type (for a full description of this 
work, see ref. |5|). 

Near the critical point the thermodynamics of the sys- 
tem (susceptibilities, correlation lengths) is determined 
by the corresponding critical exponents, which, in turn, 
are determined by the universality class of the theory. 
Thus, the universal behaviour represents a tremendous 
simplification in the effective degrees of freedom of the 
system. 

While the standard perturbative analysis can be used 
to resolve the first order nature of the transition at small 
77i ^ , it fails completely at large values of tuh — indeed, 
according to perturbation theory, the transition remains 
of first order for all Higgs masses. Thus, the physics of 
the endpoint and the universal behaviour arc inherently 
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non-perturbative. 

What kind of universal behaviour can one expect? For- 
mally, the Higgs field has SU(2) gauge x SU(2) custodia i sym- 
metry (in the SU(2) + Higgs theory, see below), but this 
remains unbroken at all temperatures. Indeed, the mass 
spectrum of the system has been investigated in detail 
both above and below the critical point, and only one 
scalar excitation (which couples to <jy<fi) becomes light in 
the neighbourhood of it. Thus, we expect the Ising-type 
universality to be realized, but also mean field-type or 
multicritical behaviour is, in principle, possible. 



II. EFFECTIVE ACTION 

An effective 3d SU(2) gauge + Higgs theory, obtained 
through dimensional reduction , accurately describes the 
static properties of the SM and many of its extensions at 
high temperatures [^) . The action of the theory is 
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The physics of the theory is fixed by the dimcnsionful 
gauge coupling g| and by the dimensionless ratios 



(2) 



where m|(/a) is the renormalizcd mass parameter in the 
MS scheme. The relations of the couplings g\,x,y to the 
full theory are computable in perturbation theory M. 
We have omitted the U(l) sector of the SM; this is justi- 
fied, since the U(l) gauge boson remains massless at any 
temperature and does not affect the transition qualita- 
tively §. 

The lattice action in standard formalism is 



S = /3 G ^(l-§Tri^) 

- H J2 3 Tr &(x)Ui(x)$(x + i) 
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FIG. 1. The phase diagrams of the SU(2)+Higgs (left) 
and the Ising (right) models. 
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Here $ is the 2x2 matrix <f> = [icr2<fr*,(f>)- The two 
actions in Eqs. (Q), (||) give the same physics in the con- 
tinuum limit a — » if the three dimensionless parameters 
0G, Ph, (3r in Eq. (||) are related to the three dimension- 
less parameters gf a, x, y in Eq. (|l|) by the following equa- 
tions iRl: 
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where £ = 3.1759115 and f = 0.08849(1). The universal 
behaviour does not depend on the lattice spacing, which 
we keep fixed through fixing flo = 5. 



III. RESOLVING THE UNIVERSALITY CLASS 

For concreteness, in what follows we shall discuss the 
universal properties of the SM by comparing it to the 
Ising model. However, one should bear in mind that the 
analysis is, by no means, limited only to verifying the 
Ising-type universal behaviour. 

The phase diagram of the 3d SU(2)+Higgs theory is 
shown in Fig. |l[ together with the Ising model phase dia- 
gram. Adopting now Ising-model terminology, let us call 
the two critical directions in Fig. [I] the Mike (perpen- 
dicular to the transition line) and the i-like (along the 
transition line) directions. Due to the lack of an exact 
order parameter, the mapping of the ft,-like and i-like di- 
rections of the SU(2)+Higgs model to the Ising model 
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FIG. 2. Top: A density plot of the 3d SU(2)+Higgs model 



at the critical point, shown on the S { 



W 2 -i) 



2 vs. Sj 



hopping 



plane. 



Each point in the plot corresponds to one measurement of the 
observables. Middle: The same as above, after a shift and a 
rotation. Bottom: A density plot of the 3d Ising model on 
the (—1 x energy) vs. magnetization plane. 



is non-trivial. This is illustrated in Fig. |2J, where the 
probability density at the critical point is plotted on the 
(Shopping, 5 f (02_ 1 )2)-plane. Only after a suitable rotation 
of the axes is the striking similarity with the Ising model 
revealed. 

This rotation closely corresponds to the rotation of the 
directions shown in the phase diagrams in Fig. The 
resulting rotated operators are conjugate to the /i-like 
and Hike couplings, and we call them the M-like ("mag- 
netization") and E-\ike ("energy") operators. Note that 
all of these operators (Shopping, S^^^ and the rotated 
M-like operator) lack the explicit M «-> —M symmetry 
present in the Ising model. This symmetry is dynami- 
cally generated at the critical point. 

However, there is no reason to restrict ourselves only to 
the two observables Shopping and Sup—xy. Any number 
of (scalar) operators can contribute to the true M-like 
and E-\ike directions. In order to improve on the projec- 
tion, it is important to consider a large number of inde- 
pendent operators. The central problem of the analysis 
is how to correctly identify the best E-\ike and M-like 
projections. 
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FIG. 3. The joint probability distributions P(Vm ,Ve) at 
the infinite volume critical point, for the volumes (a) 16 3 , (b) 
32 3 , (c) 64 3 . It is clearly seen how the distribution becomes 
more symmetric for the increasing volume, and approaches 
the Ising model distribution of Fig. ^j. 

Motivated by these considerations, we use the following 
method: 

(a) We first locate the infinite volume critical point (for 
details, see ref. [||), where all of the subsequent 
analysis is performed. 

(b) Using several volumes, we measure the fluctuation 
matrix = (stSj), s, = Si — (Si). We used up 
to 6 operators: those in Eq. (^), together with the 
operators 



FIG. 4. The probability distribution p(M, E) for the Ising 
model at the critical point (volume 58 3 ). 

X 

S L =Y,^rV^x)U t (x)V(x + i), (8) 

x,i 

where V(x) = <&(x) /\<&(x)\. 

(c) We calculate the eigenvalues \ a and eigenvectors 
V a of My. Some of the eigenvectors correspond 
to "critical" observables like M or E, and the rest 
are "trivial." They can be classified either by in- 
specting the probability distributions p(V a ) and 
p(Va,Vp), or by looking at the finite volume be- 
haviour of the eigenvalues. It is convenient to ex- 
press the eigenvalues in terms of the corresponding 
susceptibilities: Xi = -WL 3 . F° r example, the M- 
like and i?-like susceptibilities (magnetic suscepti- 
bility and heat capacity) diverge with the critical 
exponents as (L is the linear extent of the lattice) 

A M oc L^" , X E oc L a ' v . (9) 

The susceptibilities corresponding to the "trivial" 
eigenvalues remain constant. 

The method described above has much in common with 
the one used by Alonso et al M to find the E- and M- 
like directions at the critical point of the 4d U(l)+Higgs 
model, and with the method developed by Bruce and 
Wilding for the study of the liquid-gas critical point. 
Both of these rely on considering only two-dimensional 
distributions. 



IV. RESULTS 

The main part of our analysis was done on 6 lattice 
volumes, from 16 3 up to 64 3 , using the lattice spacing a = 
4/(gf/3(3) = 4/5g|. All of the simulations were performed 
close to the critical point, and the measurements were 
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FIG. 5. The divergence of \m (t°P) an d Xe {bottom) as 
a function of the lattice size. The slope of \e does not agree 
with the asym pto tic Ising value, but it is quite consistent with 
the measured Jill Ising model behaviour at these lattice sizes. 

shifted to the measured infinite volume critical point by 
histogram reweighting js). 

The fluctuation matrix analysis described in the previ- 
ous section was performed either using all of the 6 terms 
described in the item (b) above or only the 4 terms in the 
action (|3|). In both cases the projections of the eigenvec- 
tors to the original operators remained remarkably stable 
as functions of the volume. This gives us confidence that 
the method can identify the correct critical observables 
and their volume dependence. 

It is very illustrative to consider the joint probability 
distributions p(Vm,Ve)- In Fig. || these are shown for 
volumes 16 3 , 32 3 and 64 3 (using 6 operators). When the 
volume increases p(Vm, Ve) clearly approaches the Ising 
model distribution p(M,E), shown in Fig. |]. It is also 
evident that even when the SU(2)+Higgs model lacks the 
exact magnetization symmetry of the Ising model, this 
symmetry is recovered in the infinite volume limit. In 
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FIG. 6. As in Fig. ^, but for the regular observables. The 
normalization is chosen so that the average of the values is 1.0; 
the natural scale of the eigenvalues is shown in the legend. 

fact, the shapes of the distributions p(Vi, Vj) can be used 
to identify the correct i?-like and M-like eigenmodes. 

In Figs. U and || we show the behaviour of the suscep- 
tibilities Xi as functions of the volume. Only xm and 
Xe diverge as the volume is increased, the other 4 sus- 
ceptibilities do not show any critical behaviour (note the 
very restricted vertical scale in Fig. ^3] in comparison with 
Fig. ||) . The behaviour of the critical susceptibilities xm 
and Xe is compatible with the Ising model (with the 
reservation that xe does not yet attain its asymptotic 
large volume behaviour; however, it is very close to the 
Ising model behaviour at similar lattice sizes). 

The divergence of xe implies a positive value for the 
critical exponent a. This clearly excludes O(N) models 
with N > 2 (a < 0) and the mean field behaviour (a — 
0). 

The results presented above were calculated using 6 
operators in the fluctuation matrix analysis. The numer- 
ical results remain stable when only 4 operators are used, 
although in some cases small deviations from the Ising 
model behaviour begin to appear: for example, the joint 
distributions p(Vm, Ve) in Fig. || become slightly thicker 
to the E'-direction. This is in itself not surprising, since 
in the case of 4 observables the E-\ike eigenmode has the 
smallest eigenvalue, and presumably it is very sensitive 
to the quality of the projection to the operator basis f|. 
When the number of the operators is increased to 6, the 
two additional eigenvalues are smaller than the _E-like 
one. This underlines the importance of including a large 
enough number of operators in the analysis. 

In the simulations above the lattice spacing was fixed 
through Pa — Vsf a — 5. We have also located the 
critical point at /3g — 8, and Gurtler et al || have 
published results at (3q = 12 and 16. Taking into ac- 
count the 0(a)-correction calculated by Moore |1^] , we 
can extrapolate the critical coupling x c ((3g) to the con- 
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FIG. 7. The infinite volume extrapolations of x c as a func- 
tion of Pg- Giirtler et al refers to M. The O(a) -correction 
has been calculated in ref. p. 21. 



tinuum limit [3q — > oo (see Fig. ^). The final result 
x c = 0.0983(15) corresponds to the Standard Model 
Higgs mass 72(2) GeV. The effect of the U(l) gauge sec- 
tor of the SM, which was omitted here, is much smaller 
than the statistical errors quoted her e [fl3| . Since the 
experimental lower limit is ~ 88 GeV |14|, this excludes 
the existence of the SM phase transition. Nevertheless, a 
first order phase transition is still allowed in several ex- 
tensions of the Standard Model, including the Minimal 
Supersymmetric Standard Model. 



[1] K. Kajantie, M. Laine, K. Rummukainen and M. Sha- 
poshnikov, Phys. Rev. Lett. 77, 2887 (1996) ftiep-| 



ph/9605288 
[21 F. Karsch, T 



Neuhaus, A. Patkos and J . Rank, Nucl. 
Phys. B 474, 217 (1996) |iep-lat/9603004|; Nucl Phys. 
B (Proc. Suppl.) 53, 623 (1997) |hep-lat/9608087|. 



[3] M. Giirtler et al, Phys. Rev. D 56, 



(1997) § 



ep- 



Jat/970401J; E.-M. Ilgenfritz, A. Schiller and G. Strecha, 
in the proceeding s of LATTICE '98 , Nucl. Phys. B (Proc. 
Suppl.), in press [ jrep-lat/9808023[ . 



[4] Y. Aoki, Phys. Rev. D 56, 3860 (1997) ]hep-lat/9612023 
Y. Aoki, F. Csikor, Z. Fodor and A. Ukawa, in the pro- 
ceedings of LATTICE '98, Nucl. Phys. B (Proc. Suppl.) 



in press |hep-lat/9809122|; F. Csikor, Z. Fodor and J. 
Heitger, ibid. piej> ph/9809293| ]. 
[5] K. Rummukainen, M. Tsypin, K. Kajantie, M. Laine 
and M. Sha poshnikov, Nucl. Phys. B, in press [ }iep 



lat/9805013] 



[6] K. Kajantie, M. Laine, K. Rummukainen and M 
Shaposhniko v, Nucl. Phys. B 458, 90 (1996) [ aep- 
ph/9508379]; Phys. Lett. B 423, 137 (1998) [fiep- 



ph/9710538 | 
[7] K. Farakos, K 



Kajantie, K. Rummukainen, and 
M. Shaposh nikov, Nucl. Phys. B 442, 317 (1995) {iep 



lat/941209lH; M Laine, Nucl. Phys. B 451, 484 (1995) 
Hhep-lat/950400l| 1; M. Laine and A. R ajantie, Nucl. Phys. 



B 513, 471 (1998) |hep-lat/9705003] ] 
[8] K. Kajantie, M. Laine, K. Rummukainen and M 
Shaposhniko v, Nucl. Phys. B 493, 413 (1997) pep 



lat/961200q ], 

[9] J.L. Alonso et al, Nucl. Phys. B 405, 574 (1993) Jhep" 



lat/9210014] 



[10] A.D. Bruce and N.B. Wilding, Phys. Rev. Lett. 68, 193 
(1992); N.B. Wilding and A.D. Bruce, J. Phys.: Cond. 
Mat. 4, 3087 (1992); N.B Wilding, Phys. Rev. E 52, 602 



(1995) [:ond-mat/9503145[; N.B. Wi lding and M. Miiller 



J. Chem. Phys. 102, 2562 (1995) ]cond-mat/9410077|: 



N.B. Wilding J. Phys.: Cond. Mat. 9, 585 (1997) ]cond- 
mat/9610133| . 

Ill M. Hasenbusch and K. Pinn, HUB-EP-97/29 ]cond- 



mat/9706003 



[121 G.D. Moor e 



lat 



/961001E ] 



lat/970905^ ] 
[13J K 



Nucl. 
Nucl. 



Phys. B 493, 439 (1997) [ lep- 
Phys. B 523, 569 (1998) [iep- 



Kajantie, M. Laine, K. Ru mmukainen and M. 
Shaposhnikov, CERN-TH/98-281 |iep-lat/9809004| ; M. 



Laine and K. Rummukainen, in the proceedings of LAT- 
TICE '98,_N ucl. Phys. B (Proc. Suppl.), in press pep- 



lat/9809045H 



[14] The Aleph Collaboration, ALEPH 98-029, 26 March 
1998. 



5 



